Structure of the Dso{23l7) and the strong coupling 

CONSTANT QDmDK WITH THE LIGHT-CONE QCD SUM RULES 

Z. G. Wangi \ S. L. Wan^ 
^ Department of Physics, North China Electric Power University, Baoding 071003, 

P. R. China 

^ Department of Modern Physics, University of Science and Technology of China, 

O ■ Hefei 230026, P. R. China 

O . 

^ ■ Abstract 

>>: 

\ In this article, we take the point of view that the charmed scalar meson 

Dso{22>n) is the conventional cs meson and calculate the strong coupling con- 
stant QDaoDK within the framework of the light-cone QCD sum rules approach. 
The numerical values for the large scslai-DK coupling constant QDsoDK sup- 
CN ' port the hadronic dressing mechanism. Just like the scalar mesons /o(980) 

■ and ao(980), the Z)so(2317) may have small scalar cs kernel of the typical 
\ cs meson size. The strong coupling to the hadronic channels (or the virtual 
' mesons loops) may result in smaller mass than the conventional scalar cs me- 

■ son in the constituent quark models, and enrich the pure cs state with other 
. components. The 1)^0(2317) may spend part (or most part) of its lifetime as 

O I virtual DK state. 
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> 

X, 

^ , The constituent quark model provides a rather successful description of the spec- 

trum of the mesons in terms of quark-antiquark bound states, which fit into the suit- 
able multiplets reasonably well. However, the scalar mesons below 2GeV present 
a remarkable exception as the structures of those mesons are not unambiguously 
determined yet The light scalar mesons are the subject of an intense and con- 
tinuous controversy in clarifying the hadron spectroscopy, the more elusive things 
are the constituent structures of the /o(980) and ao(980) mesons with almost the 
degenerate masses. Furthermore, the discovery of the two strange-charmed mesons 
Dso(2317) and Dsi(2460) with spin-parity 0"^ and 1"*" respectively has triggered hot 
debate on their nature, under-structures and whether it is necessary to introduce the 
exotic states 0. The mass of the Dso(2317) is significantly lower than the values of 
the O"*" state mass from the quark models and lattice simulations [3] . The difficulties 
to identify the Dso(2317) and Dsi(2460) states with the conventional cs mesons are 
rather similar to those appearing in the light scalar mesons below IGeV. Those 
two states L'so(2317) and ^^1(2460) lie just below the DK and D*K threshold re- 
spectively, which are analogous to the situation that the scalar mesons ao(980) and 
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/o(980) lie just below the KK threshold and couple strongly to the nearby channels. 
The mechanism responsible for the low-mass charmed scalar meson may be the same 
as the light scalar nonet mesons, the /o(600), /o(980), ao(980) and K*(800) HEIE]. 
There have been a lot of explanations for their nature, for example, conventional cs 
states IZ1I2I1I2S1, two- meson molecular state 0, four-quark states P, etc. If we take 
the scalar mesons ao(980) and /o(980) as four quark states with the constituents of 
scalar diquark-antidiquark sub-structures, the masses of the scalar nonet mesons 
below IGeV can be naturally explained [SlEl- 

There are other possibilities beside the four-quark state explanations, for ex- 
ample, the scalar mesons ao(980), /o(980) and Dso(2317) may have bare qq and 
cs kernels in the P— wave states with strong coupling to the nearby threshold re- 
spectively, the S*— wave virtual intermediate hadronic states (or the virtual mesons 
loops) play a crucial role in the composition of those bound states (or resonances 
due to the masses below or above the thresholds). The hadronic dressing mecha- 
nism (or unitarized quark models) takes the point of view that the /o(980), ao(980) 
and Dso(2317) mesons have small qq and cs kernels of the typical qq and cs mesons 
size respectively. The strong couplings to the virtual intermediate hadronic states 
(or the virtual mesons loops) may result in smaller masses than the conventional 
scalar qq and cs mesons in the constituent quark models, enrich the pure qq and cs 
states with other component s '1 0\ ITT] . Those mesons may spend part (or most part) 
of their lifetime as virtual KK and DK states fOlE]- Despite what constituents 
they may have, we have the fact that they lie just a little below the KK and DK 
threshold respectively, the strong interactions with the KK and DK thresholds will 
significantly influence their dynamics, although the decay Dsq{2317) —>■ DK is kine- 
matically suppressed. It is interesting to investigate the possibility of the hadronic 
dressing mechanism. 

In this article, we take the point of view that the scalar mesons /o(980), ao(980) 
and Dsq{2317) are the conventional qq and cs state respectively, and calculate the 
values of the strong coupling constant QDsoDK within the framework of the light-cone 
QCD sum rules approach. The light-cone QCD sum rules approach carries out the 
operator product expansion near the light-cone ~ instead of the short distance 
a; ~ while the nonperturbative matrix elements are parameterized by the light- 
cone distribution amphtudes which classified according to their twists instead of the 
vacuum condensates O [121 • 

The article is arranged as: in Section 2, we derive the strong coupling constant 
QDsqDK within the framework of the light-cone QCD sum rules approach; in Section 
3, the numerical results and discussions; and in Section 4, conclusion. 
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2 Strong coupling constant QDsqDK with light-cone 
QCD sum rules 

In the following, we write down the definition for the strong coupling constant QDsoDK 

{K{q)D{p)\D,o{p + q)) = 9d.oDK • (1) 

We study the strong coupling constant Qd^oDK with the scalar interpolating current 
Jdso{x) and choose the two-point correlation function T^{p,q), 

T,{p,q) = tld'xe^^-{K{q)\T{j^{x)JnAO)}\0), (2) 

JdA^) = c{x)s{x), (3) 
Jj^ix) = m(x)7^75c(x) . (4) 

Here the axial- vector current Jf^{x) interpolates the pseudoscalar D meson, and the 
external K state has four momentum q with q^ = M]^. The correlation function 
Tfj.{p,q) can be decomposed as 

T^(p,g) = rp(/,(p + g)2)p^ + r,(p2,(p + g)2)g^, (5) 

due to the tensor analysis. 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach we can insert a complete series of intermediate states with 
the same quantum numbers as the current operators Josoi^) Jj^^x) into the 
correlation function T^{p,q) to obtain the hadronic representation. After isolating 
the ground state contributions from the pole terms of the Dso(2317) and D mesons, 
we get the following result. 



IAP,iP + <l))P. - {Ml-p^){Ml^^-{p + qf) 

{Ml-p^) {Ml^-ip + qr) 
where the following definitions have been used. 



+ ... 

+ ■■•, (6) 



< I J^(0) I D{p) > = if DP, . (7) 



Here we have not shown the contributions from the high resonances and continuum 
states explicitly as they are suppressed due to the double Borel transformation. 
In the ground state approximation, the tensor structure Tg (j9^, (p -|- g)^) has no 
contributions and neglected. 
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In the following, we briefly outline the operator product expansion for the corre- 
lation function T^(p, q) in perturbative QCD theory. The calculations are performed 
at the large space-like momentum regions {p + qY ^ and <^ 0, which corre- 
spond to the small light-cone distance ~ required by the validity of the operator 
product expansion approach. We write down the propagator of a massive quark in 
the external gluon field in the Fock-Schwinger gauge firstly |15j . 



(0|T{g,(xi)g,(x2)}|0)=z 

1 



-ik{xi—X'2) 



k"^ — m 



^ + m 



(2vr)4 

dvgs G'^''{vxi + (1 - f)x2 
1 



^3 



2 (P - ^2) 



A;2 



(8) 



here is the gluonic field strength, Qs denotes the strong coupling constant. 
Substituting the above c quark propagator and the corresponding K meson light- 
cone distribution amplitudes into the correlation function T^{p,q) in Eq.(2) and 
completing the integrals over the variables x and k, finally we obtain the result. 



ip + qf) - 
- Mi 



K 



du 
1 



ml — {p + uqY 



- 2 



mcg2{u) + 



Ml 

6m, 



Mu){p-q + uM^^) 



[ml — {p + ug)2]2 
1 



+ iUkMI / dv {2v - 3) / Vai^^K{a. 



{ml -[p + q{ai + vas)]'^y 



1-/3 



da- 



1 - a - /3, /3) 



{[p + (1 - Q!3 vaz)q]'^ - mlY 
1 — a — as, a^) 



dv{v — 1) 
"1 

+ I dv I das 







1 

das 

l—as 



as 



dp 



dai 







da 



(9) 



{[p+ («! + vas)q]'^ - m2}3 J ' 
In calculation, the following two-particle and three-particle K meson light-cone dis- 
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tribution amplitudes are useful fn\ IT^ [TE\ ITHl ITTj. 

< ir(g)|M(x)7^75s(0)|0 >= -zfKq^ [ du e*"^-^K(M) + x''g,{u)] 



+ fK Xf, - 



Jo 



< K{q)\u{x)tj5s{0)\0 >-- 



K 



du e™^->p(M) 



du e'"''-"(^,(M), 



< K{q)\u{x)a^^'y5s{0)\0 >= i{q^x^ - q^x^, 

< K{q)\u{x)aaf3'j5gsGf,uivx)s{0)\0 >= 

if3K[iqf^qa9uf3 - quqaQixp) " iqi,q(39ua - qvqpg,ia)] j (psKiai)e 

< K{q)\u{x)'y^'y5gsGaf3{vx)s{0)\0 >= 



iq-x{ai+va3) 



K 



qp \9a^, 

q>^ 



Xaqfj, 

q ■ X 



qa gfSf, 



q ■ X J 



■x{ai+va3) 



q ■ X J 



,iq-x{a\+va-i) 



< 



K{q)\u{x)-if,gsGap{vx)s{Q)\Q >- 



K 



qp 9a^, 



Xaq^ 

q ■ X 



qa gpf. 



+ifK-^{qaXf^-qf3Xa) / Vai(p\\{ai)e 
q ■ X J 



,iq-x(oii+va3) 



q ■ X 



(10) 



Here the operator Gai3 is the dual of the Gap, Gap = \^ap8pG^^, T^oii is defined as 
Vai = daida2da36{l — ai — a2 — a^) and 02, as) = f± + f\\ — 'f± — 'fw- The 
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twist-3 and twist-4 light-cone distribution amplitudes can be parameterized as 

<^p(m, /i) = 1 + (^30r/3 - ^p'^ C| (2m - 1) 

27 81 \ 1 



<p3K{<yi,p) 

(f)±{ai,p) 



360aia2al ^1 + A3(«i - 0:2) + 6^3^(703 - 3) 
306^{p){ai - a2)al Q + 2e(p)(l - 2a: 



Ma^,p) = 305"(p)a^(l-a3) ( ^ + 2e(p)(l-2a3; 



0||(a;j,p) = 1205^(p)e(p)(«i — a2)«i«2«3 

'1 

01l(aj,/i) = -120(5^ (/i)aia2tt3 Q + e(/^)(l - 303)^ , 
^72(m,p) = y52^p)^i(l-^i)(2^i-l), (11) 

- - - f m +m 

where Cg^, C4 and 6*2^ are Gegenbauer polynomials, e = ^1^4, 173 = "^m^'"" 
and p^ = O CSl USl EI- The parameters in the light-cone distribution 
amplitudes can be estimated from the QCD sum rules approach [T21 UHl UHl USl HI] ■ 
In this article, the energy scale is chosen to be p = IGeV. 

Now we perform the double Borel transformation with respect to the variables 
Ql = —p^ and Q\ = —{p + qY for the correlation function Tpijp', {p + q)^) in 
Eq.(6), and obtain the analytical expression for the invariant function in the hadronic 
representation, 

B,,.B,,.TpiQl,Ql) = ^<7z,.oZ,i./z,/z,.oM^,„e-*^^/^^?e-^^-.o/^^l + . . . , (12) 

here we have not shown the contributions from the high resonances and continuum 
states explicitly for simplicity. In order to match the duality regions below the 
thresholds sq and s'q for the interpolating currents J^{x) and Josoi^) respectively, 
we can express the correlation function Tp at the level of quark-gluon degrees of 
freedom into the following form, 

tAp^ (p + if) = . / . (W) 

then we perform the double Borel transformation with respect to the variables Qi = 
—p^ and Q2 = —{p+qY directly. However, the analytical expressions for the spectral 
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density p{s, s') is hard to obtain, we have to resort to some approximations. As the 
contributions from the higher twist terms are suppressed by more powers of 

-(p+q)2 ; the continuum subtractions will not affect the results remarkably, here 
we will use the expressions in Eq.(9) for the three-particle (quark-antiquark-gluon) 
twist-3 and twist-4 terms; in fact, their contributions are of minor importance. 
The dominating contributions come from the two-particle twist-3 terms involving 
the fp{u) and fa{u), we preform the same trick as Refs. [131 CHI ^^^1 expand the 
amplitudes (pp{u) and (pa{u) in terms of polynomials of 1 — 



TV 

k 



then introduce the variable s' and the spectral density is obtained. After straightfor- 
ward but cumbersome calculations, we can obtain the final expression for the double 
Borel transformed correlation function Tp{Ml,M2) at the level of quark-gluon de- 
grees of freedom. The masses of the charmed mesons are M^^^ = 2.?)VJGeV and 
Mr, = 1.865GeV^, j^J^^^^ ~ 0.45, there exists an overlapping working window for 

the two Borel parameters and M|, it's convenient to take the value Mf = M|, 
Uq = M^+M'^ ~ h ~ m'^+m'^ ' furthermore, the K meson light-cone distribution 
amplitudes are known quite well at the value = |. We can introduce the threshold 
parameter sq and make the simple replacement. 



me+ti0(l-H0)-M^ +mq (1 -up) sg 



to subtract the contributions from the high resonances and continuum states [T^ . 



^^^^^^^ -ei7-J(^M«)+ 
+e [-2fKmcg2{uo) 

+f3KMj^ / dai / ^3K{ai, 1 - ai - as, as) 2 3 

Jo Juo-ai "3 V "3 

^72 (1-%) / — da^{a,l-a- (3,(3) 

Jl^uo "3 Jo Jo 

, 2fKmMl ( f^""' da, , , , da, ['"'^'^ , n , 
-\ — / / dai / da + / dai / da 







$(a, 1 - a - as, as)]} . (15) 
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A slight different manipulation (with the techniques tak;en in the Ref.^) for the 
dominating contributions come from the terms involving the two-particle twist-3 
light-cone distribution amplitudes (Pp{u) and iPa{u) leads to the following result, 



Bm^BmiTp 



ie ' 



N 



UQ 



1— ai 



fe=0 



M 
Ml 



2\ 



1 — e A/2 



E 

i=0 



M2 



+f3KM'^ I dai / — V93i^ («!, 1 - «! - ag, as) ( 2^^^^^ — — - 3 



M2 
2fKm,Ml 



Mo J 



M2 



as 



1-/3 

/ 1 - « - /5, /5) 

'o 



/ aai / da + 

^3 J UQ—az Jo J l—UQ ^3 J UD—d'i 



dai I da 



1 — a — 03, 0:3)} . 



(16) 



In deriving the above expressions for (pp{u) + we have neglected the terms 



Mt-, here Uq 



and 



Matching the Eq.(12) with the Eqs. (15-16) below the threshold Sq and then we 
obtain two sum rules for the strong coupling constant gOsoDK, 



9DsoDK 
1 



fofDsoMoso 

d^„{u 



e "'2 ^ 



6du 

Uo 



'"c+"0(l-"0)-''-fx 

+ e 



m1+UQ(l--UQ)l\l'j, 



-2fKmcg2{uo) 



e 



+hKMl / dai I 

lO Juo-ai «3 



das . X ( ^uo - tti o 
V'3i<-(ai, 1 - - as, "3) 2 3 



03 



da-:t 



-uo a^ JQ 



as 



1-/3 



dp / (ia$(a, 1 - a - /3, /3) 



2fKm,M] 



M2 

$(a,l-a-a3, 03)]} 





dai / da + 

'^3 J UQ—Ofj, Jo J 1—uq '^3 Ju()—ag 



3 rai 
dai / 



(17) 
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9DsoDK 
1 



-e 



N 



1 _ e > ^ 

i=0 



k=0 ^ ^ 



^0 

2fKm,Mj, 



M2 
2fKm,Ml 
M2 



V93i^ (ai, 1 - ai - as, as) 2— — — - 3 



uq — ai 



as 



Q3 



dOs 
"3 



l-«o "^3 "'0 

MO 



1-/3 



dj3 / c/a<l>(a, 1 - a - /3, /3) 
das Z"^""' 



(iai / da + 



l-un "3 



«1 

(iai / da 



«0-«3 







$(a, 1 — a — as, as)} 
corresponding to the Eq.(15) and Eq.(16) respectively. 



3 Numerical results and discussions 

The parameters are taken as = (140 ± 10)Mel^, = (1.25 ± 0.10)GeK, As = 
1.6±0.4, f^K = (0.45±0.15) X IQ-^Gel/^ = -1.2±0.7, 6^ = (0.20±0.06)Ge\/2^ 

= 0.2±0.1, d2 = 0.25±0.15 O CSl 113 CHI El , Ik = O.lQOGeV, Mk = A98MeV, 
Mi5^„(2317) = 2.317GeV", Md = IMbGeV, fo = (0.23 ± 0.02)Ge1/ HSI, and /d,o = 
(0.225 ± 0.025)Ge1/ The duahty thresholds sq in Eqs.(17-18) are taken as 

So = (6.1 — 6.5)GeV^ to avoid possible contaminations from the high resonances 
and continuum states, from the Fig.l, we can see that the numerical results are not 
sensitive to the threshold parameter sq in this region. The Borel parameters are 
chosen as lOGeV^ < = M| < 20GeV^ and bGeV^ < < IQGeV^ , in those 
regions, the values of the strong coupling constant QDsoDK are rather stable from 
the sum rule in Eq.(17) with the simple subtraction, which are shown in the Figs.l- 
7. However, the values from the sum rule in Eq.(18) with the more sophisticated 
subtraction are not stable according to the variations of the Borel parameter M^. 

The uncertainties of the five parameters 5^, to'4, tus. As and 02 can not result in 
large uncertainties for the numerical values. The main uncertainties come from the 
five parameters /sx, m^, mc, fo and foaoi small variations of those parameters can 
lead to relatively large changes for the numerical values, which are shown in the 
Fig. 2, Fig. 3, Fig.4, Fig.5 and Fig.6, respectively. Taking into account all the uncer- 
tainties, finally we obtain the numerical results for the strong coupling constant, 

QD^oDK = {'^.'itlDGeV. (19) 
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\h\ 


9D,oDK{GeV) 


Reference 


0.52 ±0.17 


5.5 ±1.8 




22 




0.536 


5.68 




23 




< 0.93 


< 9.86 




24 




0.57-0.74 


6.0-7.8 




25 






10.203 




27^ 




0.88^15:- 


9.3t^:i 


This work 



Table 1: Numerical values for the parameter and the corresponding values for 
the strong coupling constant QDsoDK in the SU{?>) limit. 

The strong coupling constant Qd^oDK can be related to the parameter h in the 
heavy- light Chiral perturbation theory (201 122j , 

^——M^^-Ml\h\ 
gsp^ = ^MsMp \ 

(20) 

here the S are the scalar heavy mesons with O"*", the P are the heavy pseudoscalar 
mesons with 0~, and the vr stand for the light pseudoscalar mesons. The parameter 
h has been estimated with the hght-cone QCD sum rules |22], quark models |23], 
Adler-Weisberger type sum rule |21] , and extracted from the experimental data [23 , 
the values are listed in the Table. 1, from those values we can estimate the values of 
the corresponding strong coupling constant QDsoDK in the SU{3) limit for the light 
pseudoscalar mesons. The dimensionless effective coupling constant T/k with the 
value T/k = 0.46(9) from Lattice QCD (26j is somewhat smaller than the values 
extracted from the experimental data V /k = 0.73^*^24' here the T is the decay width 
and the k is the decay momentum. Our numerical values QDsoDK = (9.3l2j)GeV^ 
are somewhat larger comparing with the existing estimations in Refs. |22[ 123 123 
and about four times as large as the energy scale M^^o ~ 2.317GeV, and favor the 
hadronic dressing mechanism. 

Here we will take a short digression to discuss the hadronic dressing mechanism 
lOlE]; one can consult the original literatures for the details. In the conventional 
constituent quark models, the mesons are taken as quark- ant iquark bound states. 
The spectrum can be obtained by solving the corresponding Schrodinger's or Dirac's 
equations with the phenomenological potential which trying to incorporate the ob- 
served properties of the strong interactions, such as the asymptotic freedom and 
confinement. The solutions can be referred as confinement bound states or bare 
quark-antiquark states (or kernels). If we switch on the hadronic interactions be- 
tween the confinement bound states and the free ordinary two-meson states, the 
situation becomes more complex. With the increasing hadronic coupling constants, 
the contributions from the hadronic loops of the intermediate mesons become larger 
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Figure 1: The QDsoDK with the parameters and Sq . 



and the bare quark- antiquark states can be distorted greatly. There may be double 
poles or several poles in the scattering amplitudes with the same quantum number 
as the bare quark- antiquark kernels; some ones stem from the bare quark-antiquark 
kernels while the others originate from the continuum states. The strong coupling 
may enrich the bare quark-antiquark states with other components for example vir- 
tual mesons pairs and spend part (or most part) of their lifetime as virtual mesons 
pairs. 

The large values for the strong coupling constant QDsoDK obviously support the 
hadronic dressing mechanism, the Z)so(2317) (just like the scalar mesons /o(980) 
and ao(980), see Ref.^) can be taken as having small scalar cs kernel of typical 
meson size with large virtual S-wave DK cloud. In Ref . [2Z1 , the authors analyze the 
unitarized two-meson scattering amplitudes from the heavy-light Chiral Lagrangian, 
and observe that the scalar meson Dso(2317) appears as the bound state pole with 
the strong coupling constant QDsoDK = 10.203Gey. Our numerical results QDsoDK = 
{'d.?>^\)GeV are certainly reasonable and can make robust predictions. However, 
we take the point of view that the scalar meson Dso(2317) be bound state in the 
sense that it appears below the DK threshold, its constituents may be the bare cs 
state, the virtual DK pair and their mixing, rather than the DK bound state. 

4 Conclusions 

In this article, we take the point of view that the charmed scalar meson Dso(2317) 
is the conventional cs meson and calculate the strong coupling constant QDsoDK 
within the framework of the light-cone QCD sum rules approach. The numerical 
values for the scalar-Di^ coupling constant Qd^oDK are compatible with the existing 
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Figure 2: The QDsoDK with the parameters and //^-s 
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Figure 3: The QDsoDK with the parameters and 
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Figure 4: The QDsoDK with the parameters and 
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Figure 5: The QDsoDK with the parameter and fo ■ 
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Figure 7: The go^oDK with the parameter 
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estimations although somewhat larger, the large values support the hadronic dress- 
ing mechanism. Just like the scalar mesons /o(980) and ao(980), the scalar meson 
Dso{2317) may have small cs kernel of typical cs meson size. The strong coupling 
to virtual intermediate hadronic states (or the virtual mesons loops) can result in 
smaller mass than the conventional scalar cs meson in the constituent quark models, 
enrich the pure cs state with other components. The D5o(2317) may spend part (or 
most part) of its lifetime as virtual DK state. 
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